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While it is well known that neutrinos are emitted from standard core collapse protoneutron star 
supernovae, less attention has been focused on neutrinos from accretion disks. These disks occur in 
some supernovae (i.e. " collapsars" ) as well as in compact object mergers, and they emit neutrinos 
with similar properties to those from protoneutron star supernovae. These disks and their neutrinos 
play an important role in our understanding of gamma ray bursts as well as the nucleosynthesis 
they produce. We study a disk that forms in the merger of a black hole and a neutron star and 
examine the neutrino fluxes, luminosities and neutrino surfaces for the disk. We also estimate the 
number of events that would be registered in current and proposed supernova neutrino detectors if 
such an event were to occur in the Galaxy. 
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I. INTRODUCTION 

Black hole - neutron (BH-NS) star mergers are po- 
tential progenitors of short duration gamma ray bursts 
(GRBs) and have been speculated to be the site of in- 
teresting nucleosynthesis. The neutrino emission from 
the accretion disk produced in BH-NS mergers plays an 
important role in each of these scenarios. One possible 
explanation for the energetic GRBs suggests as engine a 
black hole (BH) of several solar masses accreting matter 
from a disk. Simulations of compact object mergers have 
shown the formation of such disks [TJ [5] |3J 2] . Neutrino 
transport, annihilation and losses in the disk material 
would determine the GRBs production O [7] . Several 
studies have concluded that high accretion rates would 
provide the necessary conditions for triggering GRBs 
[H [S] . Another interesting aspect is the resulting nu- 
clear products from the accretion disk around black holes 
(AD-BH) . Lattimer and Schramm [THJ E] speculated the 
possibility of a r-process in the accretion disk resulting 
of compact object mergers. Material from inner crust of 
the merging neutron star (NS), with low proton fraction, 
can be ejected from tidal tails giving place to a least a 
weak r-process. Furthermore, hot accretion disk winds 
can also produce an r-process, due to the neutrino inter- 
actions [12] . 

Given the conditions of high temperature and density 
of AD-BH, we expect a copious amount of neutrinos in 
the range of 10s of MeV to be emitted. There are several 
detectors, both in operation and proposed, that could 
register neutrinos in this energy range, including 

Super-Kamiokande (SK) [13], AMANDA [H ES], 
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KamLand [16], ICARUS [17], Ice Cube[H EDS, 
LANNDDP2], and HALO[H]. These detectors have 
been studied extensively for their ability to see a Milky 
Way supernova signal. As detection of the next Galactic 
core collapse supernova is now within reach experimen- 
tally, we can begin to speculate on future detections of 
even rarer events such as BH-NS mergers and neutron 
star - neutron star (NS-NS) mergers. Roughly, NS-NS 
mergers are three orders of magnitude more rare than 
core collapse supernovae, and BH-NS mergers perhaps 
another order of magnitude still. Given the rarity of these 
events, direct detection on a time scale of years would 
require a very large detector as we will discuss. Nev- 
ertheless it is in interesting to examine the energy and 
strength of a signal that would originate from a BH-NS 
merger and compare with the signal of neutrinos emitted 
from a proto neutron star (PNS) at the center of core 
collapse supernova. 

Previous investigations of MeV scale signals in terres- 
trial detectors of neutrinos which originate from black 
hole accretion disks have focussed on disks that might 
form from the core collapse of a rotating massive star. 
Nagataki et al, investigated the luminosity, spectrum 
and counts at SK, for neutrinos from a collapsar |20j . 
These authors used an analytical shape for the disk and 
from it derived the neutrino spectrum. Their results pre- 
dict at least one neutrino event measured in the pro- 
posed 2Mega-ton water Cherenkov detector, TITAND 
[21], originating from an AD-BH at 3 Mpc, when the 
total accretion mass, the initial mass, and the mass ac- 
cretion rate are set to 3OM , 3M Q , and O.IMqS -1 , re- 
spectively. McLaughlin and Surman have considered the 
possible distinction of neutrino spectra when neutrinos 
originated in AD-BH versus a PNS [25]. The AD-BH 
signal was found to have comparable energy spectra, but 
the disks primarily emit electron neutrinos and electron 
antineutrinos, which in addition to the timing of the sig- 
nal could produce a unique signature after neutrino flavor 
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transformation has been taken into account. 

Determining the neutrino signal from a BH-NS merger 
is more complex. In addition to calculating neutrino 
emission surfaces and estimating the effects of neutrino 
flavor transformation, general relativistic corrections are 
more important and one must determine world lines for 
the neutrinos which originate from different parts of the 
disk. In this paper we make estimates of neutrino events 
registered in several detectors for an AD-BH. We use a 
hyperaccreting disk model provided by Ruffert and Janka 
[T2l |2"5] , We also calculate neutrino luminosities and 
fluxes for the disk. We take into account general rela- 
tivity corrections for the neutrino energy and disk size. 
We discuss the consequences of neutrino oscillations in 
the fluxes and in our event counts. This paper is or- 
ganized as follows: in section [IT] we introduce the disk 



model, in section III we present the reactions included to 
calculate neutrino surfaces. In section llVl we discuss our 
results for neutrino fluxes, energies and luminosities. In 
section [V] we consider events rates in different detectors, 
in section VI we take into account neutrino mixing, and 
in section VIII we discuss our conclusions. 



formed a disk around the BH. Therefore, our analysis 
is based on one snapshot of the disk's evolution. The 
starting points for our calculation are the temperature 
T, density p m and electron fraction Y e results for each 
of the coordinates p, (f> and z of a cylindrical grid. The 
inner boundary of the model is located at p — 14 km 
and the disk surface extends until p = 300 km. While 
neutrino emission is already included in the numerical 
model, for our purposes we would like more detailed in- 
formation about the neutrinos, so in the following section 
we describe a "post-processing" of the output from this 
model. 

To compare with a steady state disk we use the model 
of Chen and Belobodorov [29]. This model is fully rela- 
tivistic. The disk is one dimensional, axially symmet- 
ric and is described by vertically averaged quantities. 
These authors worked with two different values of the 
spin parameter a — and a = 0.95, which we also use 
here. The mass of the BH is 3M Q and the accretion 
rate M = 5M©/s. The disk extension goes as far as 
p = 600km. For the vertical structure of the disk we use 
a simple hydrostatic model that assumes an equilibrium 
with the gas radiation pressure and gravity. 



II. DISK MODEL 

The results of our calculations are mainly based on a 
3D hydrodynamic simulation. However, for comparison 
purposes, we also investigate results based on a steady 
state disk. Both models are briefly described below. 

For the hydrodynamical model, we use the simula- 
tion results of a BH-NS merger by Ruffert and Janka 
[TJJ [HJ |55], for a 1.6M NS and a 2.5M Q BH with 
spin parameter a — Jc/GM 2 — 0.6 (J is the total angular 
momentum and M the rest mass of the system) . These 
authors follow the hydrodynamics of the merger with the 
Piecewise Parabolic Method [24 , including gravitational 
wave emission and neutrino emission |23j . 

In this model, general relativistic effects are included 
by using a modified Newtonian potential. The BH is 
treated as a gravitational centre surrounded by a vac- 
uum sphere. The gravitational potential $bh of the BH 
is an extension of the Paczynski-Wiita potential [3S] to 
rotating BH [27]. As function of radius r, <&bh has the 
form 
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where /3 depends on the BH spin parameter a, and rjy, 
Mbh are the event horizon and mass of the BH respec- 
tively. Artemova et al. [27] compare the results obtained 
using this potential to the results obtained with a more 
complete genereal relativistic treatment and find similar 
disk structures at 10%-20% level. 

In the model of Ruffert and Janka the Shen et al.[28] 
equation of state is used to described the stellar matter. 
The simulation was evolved until the accreted material 



III. NEUTRINO SURFACES 

Analogously to the neutrino spheres in a PNS, we can 
define the surface at which neutrinos decouple from the 
accretion disk. We follow a procedure similar to that 
outlined by Surman and McLaughlin in Ref [30J. The 
disk can be divided in regions according to their neutrino 
opacity. If in a given region the optical depth t„ > 2/3, 
then neutrinos are trapped and the disk is said to be 
optically thick. In the region where t v < 2/3 the disk 
is optically thin to neutrinos. Unlike in a PNS, neu- 
trino surfaces in a disk depend on the direction of neu- 
trino emission. We are interested in calculating neutrino 
surfaces directly above the equatorial plane of the disk. 
Therefore, we find the height of the neutrino surface, h v 
at each p, <f> by changing the lower limit in the integral 
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dz, 
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such that t v — 2/3. Here h v corresponds to a z value in 
the cylindrical grid, h max is the maximum distance in the 
z direction where matter is found and l v is the neutrino 
mean free path which is given by 
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The summation in the above equation runs over differ- 
ent neutrino scattering process, which we describe below, 
n k — n k(p, 07 z) is the associated particle density of each 
process and 
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is the corresponding cross section averaged over the 
Fermi-Dirac flux 



<f>(E u ) 
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with g v = 1 and with an assumed neutrino chemical po- 
tential /i„ = 0. Strictly speaking the flux (j>(E v ) is also a 
function of T and therefore a function of p as well as of 
and z. We will point out this dependency when needed 
by noting the value of T at which <f> is evaluated. 

As matter is dragged into the black hole the medium 
becomes hotter and denser, and nuclei dissociate. There- 
fore, we consider neutrino scattering from protons, neu- 
trons and electrons. We have the current charged reac- 
tions for electron (anti)neutrino (z7 e )^ e : 
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v e + P — > e. + n, 
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and for all (anti) neutrino flavors the neutral current pro- 
cesses 



v + p — > v + p, 



v + n — > v + n. 



v + e 
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The cross section for neutrino absorption (Eq, [6| , includ- 
ing weak magnetism effects Wm [31] is 
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where 
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gA = 1.93, m^fn) is the electron (neutron) mass, A = 
1.23 is the neutron proton mass difference, and cto = 
AG 2 F m 2 / nh A . For scattering of electron antineutrinos 
from protons, Eq. [7] we have 
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with 



W M = 1-7.1 



(15) 



The cross section for the neutral current process of Eq. 
[8] valid for all neutrino flavors is [32] 
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with Cy = 1/2+2 sin 2 Qw, the Weinberg angle sin 6^ = 
0.23 and Ca — 1/2, whereas the corresponding expres- 
sion for Eq[9]is 
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The same expressions hold for antineutrinos if we change 
gA = —gA- This change does not have consequences in 
the above expressions. 

In order to obtain a(E v ) for the electron neutrino elas- 
tic scattering from electrons (Eq. 10 1, we integrate the 
differential cross section 133 
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where, 
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and evaluate it between the maximum and minimum 
values of T e , corresponding to T e = m 2J ^ 2 _E and 0, re- 
spectively. For v x (x — [J., t) scattering from electrons 
we have the same expression with Ca = — 1/2 and 
Cv = 2 sin 2 <dw — 1/2. Antineutrino scattering requires 
changing C A = -C A - 

Finally, an approximate form for the cross section for 
neutrino-neutrino annihilation, Eq |ll[ is taken from Di- 
cus [M] and Goodman [35], assuming E V ,E^ » m e [29] 
(for the complete expression see [36"]). 
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for antineutrinos. Here (E u t v \) are the average 
(anti) neutrino energies and 
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To calculate the proton n p and neutron n n number 
densities needed in Eq. [3] we start from the electron 
fraction fraction Y e and the mass density p m of the disk 
model. We assume charge neutrality Y e — Y p so we have 
rip = p m N A Y e and n n = p m N A {\ -Y e ), where N A is 
the Avogadro's number. We find the electron number 
density n e using 



p e - + p e + = 0, 



(24) 



and Fermi-Dirac distributions for electrons and positrons 
with the corresponding grid values of temperature. We 
proceed by finding /x e - such that n p = n e ~ — n e +, and 
with this value of /i e - we find from the Fermi-Dirac dis- 
tribution that characterizes n e for each grid point. 

The procedure described above allows us to find the 
temperatures T v at which neutrino decouple. Those cor- 
respond to the temperatures at h v in Eq. [2] We show 
the result of this calculation in Fig. [I] This figure shows 
a transversal cut of the resulting neutrino surface corre- 
sponding to a polar angle <p = 20° in the original numer- 
ical simulation. 

The first fact to be noticed is its irregular shape, in con- 
trast to the symmetric one of a neutrino sphere. Also, 
because the material in the disk is relatively neutron rich, 
v e absorption on neutrons has a more significant contri- 
bution to the mean free path, than v e absoprtion on pro- 
tons. Thus the electron neutrinos decouple at the lowest 
temperatures. The mu and tau type neutrinos and an- 
tineutrinos lack these charged current interactions and 
decouple at the highest temperatures. This is the same 
hierarchy of energies that is seen in neutrinos emitted 
from the PNS. 

Fig. [2] shows the electron antineutrino surface temper- 
atures for the whole disk. The disk is in the x, y plane 
and the color scale represent the temperatures. Blue 
corresponds to T = 0. The hotter D e are closer to the 
BH (black circle in the center). These temperatures are 
higher than the PNS. While the core temperatures in 
a PNS arc of the order of hundred MeV, the tempera- 
tures of the AD could be around 20 MeV. However, the 
density profile of the disk is very different from that of a 
PNS. When a neutrino decouples in a PNS it has diffused 
through a denser medium which is assumed to be more or 
less symmetric in all directions, giving as a result a spher- 
ical shell for the neutrino surface with lower temperature. 
In the case of the disk, as matter is dragged to the BH, 
a funnel is formed around the BH vicinity. Therefore, 
changes in density are more abrupt. The medium is less 
dense close the BH, increases rapidly between 30 to 60 km 
and then decreases again as r increases. Neutrinos emit- 
ted close to the BH travel through a less dense medium 
for less time, resulting in higher temperatures compared 
to those in a PNS. 

Fig. [3] shows a 3D image of the D e surface. The height 
represents h u while the color scale represents T Pc . The 
biggest contribution to the antineutrino flux comes from 
regions closer to the boundary with the BH. There, the 
h v is smaller which translates in higher decoupling tem- 



peratures. As described in the text, this is a consequence 
of the density profile. Regions where h v is high and far 
away from the center contribute less to the flux. 
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FIG. 1: (Color on line) Transversal cut of neutrino surfaces 
at 4> = 20°. The solid line corresponds to electron neutrino 
surface whereas the dashed and dotted lines correspond to 
electron antineutrino and tau neutrino respectively. The circle 
around r = represents the BH boundary. 



IV. FLUX, LUMINOSITY AND ENERGY 

Using the values of T v we calculate fluxes, luminosities 
and averaged energies for each neutrino flavor. For the 
neutrino luminosity we integrate E v over the neutrino 
surface, assuming the Fermi-Dirac distribution of Eq. [5] 







L v = / d<; 


b ( pdp ( 


Jo 


J Pmin JO 



E v <l>(E v ,T v )dE v . (25) 



Here p m in = 14 km, p max corresponds to the boundary, 
in the radial direction, on the x, y plane, of the opti- 
cally thick region, and <f> is the polar angle. L v is the 
total energy emission rate and does not depend on a spe- 
cific observation point. The total luminosity is obtained 
by adding the resulting luminosity for each neutrino fla- 
vor. Similarly, to estimate the total number of neutrinos 
emitted per unit time /, we integrate the flux over the 
neutrino surface, 
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Our estimate for the average neutrino energy is then 
given by 
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FIG. 2: (Color on line) Temperature profile of the electron 
antineutrino surface seen from the z axis. The black frame, 
x = [—70,75], y — [—75,70] km, encloses the antineutrino 
surface. The temperature scale (on the left) goes from blue 
T — 2 MeV, to red T ~ 14 MeV. The black circular area 
represents the black hole boundary, r = 2r a . 




T [MeV] 
15.00, 



FIG. 3: (Color on line) Electron antineutrino surface seen at 
some inclination angle (see the x, y, z axis on the lower left 
corner). The height corresponds to h v as in Eq. [2] The color 
scale corresponds to the neutrino temperatures, also shown in 
Fig. [2] The black area in the center represents the boundary 
with the BH, r = 2r s . 



Because of the presence of the black hole the above 
quantities will differ from those measured by an observer 
at infinity. A source located at a distance r from a black 
hole emitting at energy E will be observed at infinity to 
have an energy, [55] , 



E* = 



E 
1 + z' 



(28) 



The redshift factor 1 + z consists of a Doppler part and 
a gravitational part (see Ref . [37] for a derivation) . The 
Doppler term depends on the ratio of f2, the angular 
velocity of the emitting gas, to the speed of light. We 
find that the Doppler term is several orders of magnitude 
smaller than the gravitational term. Therefore, we use 
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and then the energy observed at infinity is 
E* = \g 00 \ 1/2 E, 



(29) 



(30) 



where goo is determined by the space time metric. In the 
case of a non-charged, rotating black hole the curvature 
of the line element can be written in the Kerr geometry 
as ESI, 



ds 1 



-(A/£ 2 ) [di-a 2 sin 5 



(sin 2 9/i 2 ) [{r 2 + a 2 )dct> - adt] 2 



(£ 2 /A) dr 2 + fete 2 
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where 

A = r 2 — r s r + a 2 

£ 2 = r 2 + a 2 cos 2 6, (32) 
and the Schwarzchild radius r s = 2M . goo is given by 
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Note that the parameter a given above is the spin of the 
black hole, not the disk. Before the merger the spin pa- 
rameter of the black hole is a = 0.6. The Kerr metric re- 
duces to the Schwarzschild case when a = 0. We will use 
the Kerr metric to describe the neutrino energy redshift. 
However, for simplicity (and because the spin parameter 
introduces only a small correction unless it is nearly one) 
we use the Schwarzschild metric to calculate corrections 
due to the neutrino ray bending. An observer at infinity, 
from a BH will detect a luminosity L* — \goo\L from an 
object which has a luminosity L [39] . For our disk we 
calculate the redshift in energy |3oo( r ^)l as i n Eq 33 at 

„2 



the point of decoupling r 2 = h 2 + p 2 , with h v and p v 
the corresponding cylindrical emission coordinates. Here 
again L* is the total emission rate and does not take into 



account a specific location of the observer 
have, 



.Then we 



6 



TABLE I: Observed E*(L*) and emitted E(L) averaged neu- 
trino energies (luminosities) for an AD-BH and for a PNS. 

E{McV) E*{McV) E(MeV) L(ergs/s) L* (ergs/s) 
Disk Disk PNSgl] Disk(xl0 53 ) Disk(xl0 53 ) 
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for the observed flux, luminosity and averaged energy 
respectively. Table [I] shows our results for averaged en- 
ergies and luminosities, both emitted and observed. For 
comparison we have added the values corresponding to 
a PNS. The total observed luminosity for this AD-BH is 
1.6 x 10 54 ergs/s, a hundred times larger than the lumi- 
nosity of a PNS, L ~ 10 52 ergs/s. 

We find the results for the Chen and Beloborodov 
steady state disk, applying the same technique for the 
calculation of neutrino surfaces, energies and luminosi- 
ties. The energies obtained for spin parameters, a = 
and a = 0.95, are lower than the results of Table [T] For 
example, the average energy for electron antineutrino 
is E v = 11 MeV and E v = 14.4 MeV, for a = and 
a = 0.95 respectively. This difference can be understood 
in terms of the temperature and density dependence with 
the distance to the BH. For both models the tempera- 
tures and densities are similar in the region close to the 
BH, where most neutrinos are emitted. The Janka and 
Ruffert model predicts highly fluctuating temperatures 
and densities that drop faster with distance. The Chen 
and Belobodorov model leads to symmetric neutrino sur- 
faces and a smooth temperature decay. Therefore, as the 
distance increases neutrinos with low temperatures con- 
stitute a significant fraction of the final spectra lowering 
the average energies. 



V. NEUTRINO SPECTRA 

So far we have determined neutrino energies and lumi- 
nosities, both emitted and observed at infinity. In this 
section we obtain the neutrino fluxes as seen at a fixed 



point and apply our results to estimate the number of 
neutrinos registered at that specific point. 

The number of neutrinos emitted per unit energy, per 
unit area, per second, reaching an observer located above 
the disk plane, on the z axis, at a distance z eva is given 

by 
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Here 9 ma x is the maximum angle formed by the outer 
edge of the neutrino surface and the z axis, as seen by an 
observer at z eva . E* is the redshifted energy. 

If the BH did not affect the neutrino trajectories the 
angle max subtended by an observer at z eva would be 
given by ttm(0 max ) = r max /z eva , with r max defined by 
the boundary of the neutrino surface. However, the pres- 
ence of the BH does bend the trajectories. To take this 
into account, we can follow backwards the neutrino tra- 
jectories by tracing null geodesies leaving the observer 
at z eva and reaching the disk [42 . The angle subtended 
by the disk according to the observer can be calculated 
in terms of the impact parameter b, which is a con- 
stant over the trajectory. The diagram in Fig. |V| shows 
the effect of the BH on the neutrino trajectories. Neu- 
trinos leave the neutrino surface at the emission point 
r„ = (h 2 v + pi) 1 / 2 - Their trajectories bend according 
to their separation from the BH. The influence of the 
gravitational field is less strong when the neutrinos are 
far from the BH. b can be visualized by assuming that 
when the neutrino is far away from the BH it travels in 
a straight line. The impact parameter is the distance 
between the closest approach of the continuation of this 
straight line and the center of the BH [43] . 

According to the definition of b we have, 



sin 9 



(38) 



and therefore for the element of solid angle seen by an 
observer at infinity d£l — bdbdtfi/ z\ va [U]. On the other 
hand, b that a neutrino going straight up will encounter 
at the emission point r v is in the Schwarzschild metric 
given by 



b = 



(l-r./r„)V2 



and 



db- 



dr„ 



(1 - r s /rv)V2 



2r(l-r,/r„) 



(39) 



(40) 



Using this expression, the integral over b can be written 
in terms of r v (or equivalent p v ) and then the effective 
flux will be written as 



4tt J{i 



bdb 



d(b x <p(e:). (4i) 



with limits given by boundaries of the neutrino surface. 
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FIG. 4: Diagram showing the effects of the BH gravitational 
field on neutrino trajectories. The black dot represents the 
BH, r v is the emission point on the neutrino surface, which 
we have sketched with the horizontal line, b is the impact 
parameter of the trajectory (thick line) and 9 is the angle 
measured by an observer at z eva . According to this observer 
the emission point would be located at r*. 



A. Neutrino Counts 

We now use the neutrino fluxes calculated as above to 
estimate the number of counts in current and proposed 
neutrino detectors. We assume that we observe the disk 
from the z axis at distance z eva = 500 km and then we 
re-scale our results to lOkpc. We evaluate the count rate 
R, in a given detector, by integrating the effective flux 
over the energies recorded at the detector 



R = N T 



4> en (E u / 9m (r))a(E*)dE* 



(42) 



Here Nt is the number of targets in the detector 
and E* and E are related by Eq. [28] The flux 
(f> e ff( Ev/ gQo(r)) is the redshifted neutrino distribution 
of Eq|4l] with 4>{E y ) evaluated at T„, E* is the detected 
energy, and E t h and cr(E*) correspond to the threshold 
energy and cross section of the detector. Using Eq. 41 
we write R as 



R = 
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eva J Et 



cj>(E„/g 00 (r))a(E*)dE* 



(43) 

The cross sections in Eq. [43] depend on the detector 
under consideration. For SK this corresponds to the ex- 



pression Eqjl4| The energy threshold in SK is 5 Mev. 
Our results are for a 32 kton volume. 

We can easily re-scale our results to the bigger volumes 
of the proposed detectors UNO (580 kton) g51 HB] and 
Hyper-K (1 Mton)|37]. On the other hand, the large scale 
AMANDA (in operation) and Ice-Cube (in construction) 
detectors, though designed to detect high energy neutri- 
nos, have been discussed to detect supernova neutrinos 
if the number of them is large enough to allow extract- 
ing their signal from the detector background [HI ITS] . 
We think this is the case for the AD-BH neutrinos and 
we estimate counts in these detectors using E t h — and 
an effective volume of 414 m 3 for each of its 680 optical 
modules (OM). Ice-Cube is an extension of AMANDA 
with 4800 OM. 

For supernova neutrino detection, Argon based facili- 
ties such as ICARUS or LANNDD rely on the reactions: 

Ar^e-+ 40 K*, (44) 



40 



and 



40 



Ar^e++ 40 Cl*, (45) 

as well as elastic scattering from electrons. We estimate 
our counts in these detectors assuming a fiducial volume 
of 3 kt and 70 kt for ICARUS and LANNDD, respec- 
tively. We take the cross sections from Refs. [1H1 US] • 

In a lead based detector, such as the proposed HALO, 
the high number of neutrons Pauli-block the electron an- 
tineutrino absorption. However, the neutrino scattering 
cross section per nucleon is large and most of the events 
produce neutrons. Neutral current process also occur 
with the emission of neutrons. We have used the cross 
sections by J. Engel et al. [50 to calculate rate counts 
for the charged current reactions: 



v.. 



208 



Pb 



207 



208 



U e -K" J Pb -^ 20G Bi + 2n 
and for the neutral current processes, 



v 



208 



and 



Bi + n + e~, (46) 

+ e~, (47) 

t processes, 

Pb -^ 206 Pb + 2n, (48) 

v+ 2QS Pb^ 207 Pb + n. (49) 

KamLAND also offers the possibility of neutrino detec- 
tion from supernovae via elastic scattering from protons. 
Beacom et al. |16j predict neutrino proton recoil spec- 
trum for those processes including quenching of the pro- 
ton scintillation light and different detector backgrounds. 
We follow their treatment for the proton energy thresh- 
old. We start from the differential neutrino-proton cross 
section 



da 



dT n 



8m2 



(C v + CaY + (cv - c a y 1 - 



T, 



<c 2 



c 2 



TflpTp 



(50) 
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where m p and T p are the mass and recoil kinetic energy of 
the proton respectively, Cy = 1/2 — 2 sin 2 Q\y and Ca = 
1.27/2. The same expression holds for antineutrinos with 
Ca — ~Ca- For a neutrino of energy E v the maximum 
proton kinetic energy T™ ax is 



rpmax 

P ~ 



2Ef, 



2E V 



(51) 



In a scintillator the light output from low energy protons 
is reduced relative to the light output for an electron de- 
positing the same amount of energy. Taking into account 
this proton quenching the threshold on the proton kinetic 
energy in KamLAND is T™ m =1.2 MeV [IB]. Therefore 
we integrate the differential cross section as 



T p m< " da, 



Up 



fmin dTp 



(52) 



and when replaced in Eq[43] a is integrated with respect 
to E v with 



Eth 



m p T™ ln 



-1 1/2 



(53) 



We obtain the total counts by multiplying the rate 
count R by the duration of the signal. This time is de- 
termined by the amount of the total binding energy, Eb, 
that can be emitted by neutrinos. We note this quantity 
as E B . McLaughlin and Surman [22\ estimated that from 
the total binding energy at r s , E B = 9(M/M ) x 10 53 
ergs (M being the NS mass), 20% is released in the 
form of neutrinos. This means E B ~ O.lMc 2 . Other 
estimates of accretion onto BHs result in a budget of 
GM/r ms ~ 0.1c 2 = 10 20 ergs/g, with r ms the radius of 
the marginally stable orbit |51j . or, if a relativistic disk 
accretion onto a Schwarzschild BH is used 5.7% of the 
rest mass energy [3"9"| . 

Roughly speaking, neutrinos would be emitted during 
an interval of time dt = E B /L, with L the total neutrino 
luminosity. However, due to the strong gravitational field 
caused by the BH, observers on Earth will detect a longer 
signal lasting an interval dt* , 



dt* = 



dt 



l5oo I 1 / 2 



(54) 



To estimate an overall l^ool 1 ^ 2 , we use our results for 
the observed and emitted energy averages and the re- 
lation between emitted and observed energies, Eq. [30] 
Using the estimates for E B , we find that observers on 
Earth will detect a signal lasting for an interval of time 
dt* = 0.15 — 0.07 s. The estimate of the neutrino sig- 
nal for a supernova is 10 s and for AD-BH from NS-NS 
mergers ~ 1 s [55]. Therefore the AD will emit a signal 
approximately 100 times shorter but, according to the 
results in table [I] more luminous. These fact offers the 
possibility of distinguishing the AD-BH spectrum from 
that of PNS. 



TABLE II: Neutrino counts from a BH-NS merger as regis- 
tered at several facilities. Lines between rows separate de- 
tectors according their principle of detection. The size of the 
detectors are indicated in parenthesis. We include the PNS 
count estimates for some detectors also indicated in parenthe- 
sis. 





Ve + p — > n + e T 


v + e — > v + e 




SK(32 kton) 


9100 


390 




UNO(580 kton) 


165000 


7100 




T T TV / 1 TV if J- \ 

Hyper-K(lMton) 


no d r\r\r\ 

284000 


12280 




Amanda(680 OM) 


74000 


2800 




iceOube(4oUU UMJ 


O2250U 


202U0 




tr IN b(orv) 


8300 


32U 






V + p — > V +p 






KamLAlNL) (1 kton) 








PNS 


273 








^ e + 208 Pb 


^+ 208 Pb 


total 




^207(6) Bi+e 


_ > 207(6)p b 




HALO (80 ton) 


24 


23 


47 


PNS 






43 




^ e )+ 40 Ar 


v + e — > v + e 






- e(e+) 






ICARUS (3 kton) 


331 


30 




LANNDD(70 kton) 


7700 


700 




PNS (ICARUS) 


203 


41 





Table ILT1 shows our results for the 

BH-NS neutrino counts registered at the several detec- 
tors discussed in this section, with an estimated signal 
time dt* = 0.15s. Unless indicated we have added results 
over all neutrino flavors, and/or over the final products. 
For example reactions on ICARUS correspond to the sum 
of counts for absorption of v e plus v e . For comparison 
purposes we also show the current estimates for a PNS. 
From table [H] we can see that, in most of the cases, an 
AD-BH located a 10 kpc will generate more counts that 
a PNS at the same distance. This is general for all the 
detectors considered here. In particular SK will record 
1000 more events from 

a BH-NS merger. 

Using the neutrino surfaces of the steady state disk 
model we estimate the counts obtained in SK. We pre- 
dict an event rate of 8000/s and 31000/s for a = and 
a = 0.95 respectively. The higher rates are a result of 
the higher temperatures of the rotating disk. To esti- 
mate the total counts at a given detector we need to 
estimate the binding energy available for neutrino emis- 
sion. As in the previous case, to calculate the duration of 
the signal, we assume that the total mass dragged to the 
BH corresponds to 1.6Mq and then divide by the total 
neutrino luminosity. This give us a total of 5400 counts 
in SK for a = and 4000 for a = 0.95. The total num- 
ber of counts in the rotating disk are less because such 
disk is more luminous and therefore the available binding 



9 



energy radiated as neutrinos would be consumed faster 
than in the case of a non-rotating disk. The extent to 
which steady state and dynamical models can converge 
on similar structure is not yet clear. We continue our cal- 
culations with the multidimensional model. As the disks 
that occur in compact object mergers become better un- 
derstood, the methodology presented here can be applied 
to these new disk models as well. 

Re-scaling to different distance we can use these re- 
sults to determine how far we can see neutrinos from 
NS-BH accretion disks. For example, to detect neutrino 
from a supernova in AMANDA a vast amount of counts 
is needed, so the signal stands out from the background 
noise [14] . If we speculate that for being able to distin- 
guish the signal from BH-NS merger, then it is necessary 
at least the same amount of counts as from supernova at 
10 kpc, then AMANDA could reach AD-BH as far as 30 
kpc. This estimate assumes a window interval of 10 sec; 
however the BH-NS signal lasts around 0.1 sec, therefore 
this 

is a conservative estimate. A more optimistic figure is 
obtained if we consider a large scale detector like UNO, 
which would detect at least one count from AD-BHs lo- 
cated as far as 4 Mpc. Compare this result with Advance 
LIGO, whose reach for NS-NS is 300 Mpc and for NS-BH 
is 650 Mpc [53]. In order to see a NS-BH merger in the 
same distance range, but in neutrinos, requires a Gigaton 
scale detector. 



VI. NEUTRINO MIXING 

In their way from the accretion disk toward a detector 
on Earth neutrinos will go through flavor transforma- 



U = 



The coefficients of U depend on the mixing angles 
#12, #13, and 6*23, and the CP- violating phase 8. In U, 
Cij = cos 9ij and Sij = sin 0y . So far there is not a way 
to discriminate the neutrino mass ordering of the neutri- 
nos. The ml < m\ < m 3 relation is known as the normal 
hierarchy (NH) while m§ < m\ < m| is refered to as the 
inverted hierarchy (IH) . We consider oscillation scenarios 
for both NH and IH. 

For a supernova density profile there is one matter res- 
onance at low densities, usually called the L resonance, 
and another matter resonance for high densities, the H 
resonance. The L resonance mixes the mass states v\ and 
u 2 involving the mass difference 5m\ 2 , which can be de- 



tion. At present there is some uncertainty in the flavor 
transformation, in part due to unknown neutrino param- 
eters such as the hierarchy and the third mixing angle, 
and in part due to the lack of a complete calculation of 
neutrino flavor transformation for neutrinos leaving ac- 
cretion disks. However, since most detectors measure fla- 
vor dependent signals, the neutrino oscillations will have 
some effect on the number of counts. In this section we 
briefly review the possibilities for flavor transformation, 
and then take a few scenarios to demonstrate the range 
of possibilities. We follow the procedure of Kncllcr et al. 

Neutrinos have mass and therefore the weak states 
e, /i, t can be described as linear combinations of the mass 
states mi, ni2, 17*3. In the presence of matter, the Hamil- 
tonian describing the evolution of the neutrino states is 
neither diagonal in the flavor basis nor in the mass ba- 
sis. The coefficients describing the linear combination in 
any basis will oscillate, and their behavior will depend 
on the medium through which the neutrinos propagate, 
their energy, the differences between the squares of the 
masses <5m 2 „- = mf — m 2 and the mixing angles conncct- 
ing the flavor and mass basis. The matrix U relating the 
flavor and mass states can be written as 



(55) 
I 

termined by solar neutrino experiments <5m 2 2 ^ Sitiq and 
9i2 ^ 9q. Because of the uncertainty in the mass hier- 
archy the H resonance can involve different mass states. 
In the normal hierarchy the states mixed are v 2 and is 3 , 
and involves the mass splitting <5m 2 3 which is unknown. 
In the inverted hierarchy the mixed states are v\ and v 2 
involving the mass splitting (5m 2 3 . For both hierarchies 
the mixing angle is the small #13. The experimental limit 
to 

date is sin 2 29\ 3 < 0.19 [55]. The probability for a 
neutrino going from one matter state to the other in L 
resonance is denoted here as P^(E V ), and in the H res- 
onance as Ph{E v ) and Ph(Ep). Antineutrinos do not 



C12C13 S12C13 si 3 e~ %s 

-S12C23 - Ci 2 Si 2 S 2 3e lS C12C23 - Sl2Sl3S23e l<5 C13S23 

\ S12S23 - c 12 s 13 c 23 e lS -c 12 s 23 - s 12 s 13 c 23 e l5 c 13 c 23 
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cross from one state to another in the L resonance. If 
the density changes slowly the neutrinos can propagate 
adiabatically. Therefore a neutrino produced in a 

matter state will remain in the same 

matter state as propagates through the medium. This 
means the crossing probability is close to zero and the 
dominant flavor state after passing through the resonance 
would have changed. The resonance is then "adiabatic" . 
A "non-adiabatic" resonance corresponds to a crossing 
probability closer to one. 

Neutrino self-interactions can also affect the probabil- 
ity of detecting one flavor neutrino or another, and this 
is important in the region relatively close to the emis- 
sion point where these self-interactions are large. For a 



review see [56; . As described in the review, depending 
on the density profile, the luminosities of the neutrinos, 
and the neutrino parameters, e.g. the hierarchy, these 
self-interactions can produce a range of new behaviors. 
These include a spectral swap, a spectral split, and no 
change at all. We follow the expressions in Ref . [54] , and 
denote the survival probablities from this region where 
self-interactions can dominate as Psi(E v ) and Psi(Ep). 

While ideally one should follow the evolution of the 
neutrino wave functions fully as they travel out of the 
accretion disk, as long as the relevant, H, L and SI re- 
gions are well separated we can approximate the survival 
neutrino and antincutrino probabilities, p, p, as follows. 

In the normal hierarchy, 



P = [\U el \ 2 P L + \U e2 \ 2 (l - P L j\ [P H (1 - Psi) + (1 - Ph)Psi] + \U e3 \ 2 [PhPsi + (1 - Ph)(1 ~ Psi)] 

P = \U el \ 2 [PhPsi + (1 - Ph){± - Psi)] + \U e z\ 2 [P H (l - Psi) + (1 - Ph)P S i] , (56) 

and in the inverted hierarchy, 

V = [\U el \ 2 P L + \U e2 \ 2 (l - P L )\ [P H P SI + (1 - P H )(1 - Psi)} + \U e3 \ 2 [P H {\ - P SI ) + (1 - P H )(1 - P SI )} 

V = \U el \ 2 [P H (1 - P SI ) + (1 - P H )Psi] + \U e3 \ 2 [PhPsi + (1 - Ph)(1 - Psi)] ■ (57) 



The observed fluxes, after flavor oscillations, are ex- 
pressed in terms of the survival probabilities as 





Ki-p)<.. 




(58) 




Ki-p)4>L 




(59) 




-p)< + (i- 




(60) 


> 






(61) 



where <fi° are the initially produced fluxes and x is cither 
of the heavier flavors, x = r = /i. Similarly, using the 
initially emitted energies the results of flavor oscillations 
on the energies. We insert the corresponding fluxes, in- 
cluding oscillations, in Eq. 



43 and calculate the counts 



at the detector. We consider oscillations scenarios with 
both normal and inverted hierarchies. For simplicity we 
calculate the survival probabilities assuming completely 
adiabatic or non-adiabatic resonances, that is Pl = or 
Pl = 1, and similarly for Ph and Psi- Combining all 
the values of Ph and Pl we have a total of 8 scenar- 
ios. Finally, we use sin 2 612 = 0.311 and sin 2 #13 = 10~ 4 
|55j for the mixing angles needed to calculate the matrix 
coefficients EZy. 

From the different scenarios studied we chose the ones 
which give extreme values. Therefore, we can have an 
idea of what could be the upper and lower limit of our es- 
timates and how much oscillations can affect them. Sev- 
eral scenarios reproduce the same values. We pick two 
for our discussion, SI and S2. SI correspond to a com- 



pletely adiabatic H resonance (Ph — Ph — 0) and non- 
adiabatic L resonance (Pl = 1) in the NH. S2 considers 
IH, with P L = Ph = Ph = 0. We have used P SI = 0, 
and Psi = 1 for both NH and IH. Table |lHl shows the 
results of neutrino oscillations for energies and number 
neutrinos emitted per unit time /. Table IV shows the 



effects of neutrino mixing in the registered events for the 
same oscillation scenarios. The energies of v e and v e 
are larger while energies for v x and v x decrease. As a 
consequence of the neutrino mixing the energies of all 
neutrino flavors get closer together and the values of the 
fluxes also become closer together. Counts in most de- 
tectors are larger as a result of neutrino mixing, except in 
KamLAND where our results remain constant. The extra 
counts in detectors where initially an electron neutrino or 
antineutrino is scattered can be easily understood: there 
is an extra contribution to the flux from (v x )v x emitted 
at high temperatures from the neutrino surface, see table 
III 



VII. SUMMARY AND CONCLUSIONS 

The physics of NS-BH mergers is very important for 
our understanding of GRBs and the production of rare 
isotopes. 

Using a state-of-the-art model for a NS-BH merger 
we have calculated neutrino surfaces, fluxes luminosities 
and energy averages for the neutrinos emitted from the 
emerging AD. Neutrino surfaces as seen from the z axis 
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TABLE III: Effects of neutrino mixing in the number of neu- 
trinos emitted per second / and in energy. We have included 
the values before mixing for comparison (NOsc). 

NOsc SI S2 



^IVlC V J 








E Ve 


17.3 


20 


26 




23.4 


24 


26 




26 


25 


22 




26 


25.6 


25 


(xl0 57 ^/s) 








Jv e 


6.0 


6.4 


7.1 


fv c 


6.5 


6.7 


7.1 


.fl'.r 


7.1 


6.9 


6.6 




7.1 


7.0 


6.8 



TABLE IV: Effects of neutrino mixing in events registered at 
several facilities. We have included the values before mixing 
for comparison (NOsc). For ICARUS we only show counts 
due to the current charged channel (CC) 





NOsc 


SI 


S2 


v e + p — > n + e 4 " 








SK 


9100 


9800 


11460 


AMANDA 


74000 


79970 


93200 


v + e — > v + e 








SK 


390 


430 


490 


ICARUS (CC) 


331 


450 


710 


HALO 


47 


62 


91 


KamLAND 


470 


468 


467 



have an irregular shape. Due to large cross sections the 
surface for v e extends higher above the disk plane com- 
pared to the other neutrino flavor surfaces. We find that 
the redshifted neutrino energies are higher than in the 
case of a PNS. Our estimates are E Ue p=> 17.3 < E^ e ps 
23.4 < E Vx ps 26 MeV, which include redshift due to "the 
strong BH gravitational field. The AD-BH neutrino lu- 
minosity from the NS-BH merger is also larger than the 
one of a PNS. We find a total neutrino luminosity ~ lfj 54 
ergs/s; this is two orders of magnitude larger than the 
luminosity of a PNS. According to previous works this 
luminosity could be enough to trigger short GRB [5]. 

These results are based on a snapshot of the disk evolu- 
tion. We have calculated the same quantities for a steady 
state model and found that the average energies are lower 
compared to the hydrodynamical model. This is due to 
the contribution of low temperature neutrinos emitted at 
large radial distances, p, from the disk to the spectra in 
the steady state cases. We explored the influence of the 
spin parameter in the steady state disk using two differ- 
ent values. We have found that higher spin parameters 
lead to higher neutrino energies and luminosities. As a 
consequence the neutrino counts per second are larger 



for the rapid rotating disk. However, if the fraction of 
binding energy released in neutrinos were the same for 
both, rotating and non-rotating disks, the signal in the 
rotating disk will be shorter resulting in less total counts. 
The order of magnitude of the luminosity for the steady 
state disk with a=0.95 is the same as for the 3D model 
and as it has been discussed could be a good candidate 
for triggering GRB [5]. 

As part of our estimates, we have explore the conse- 
quences of neutrino mixing on the neutrino energies and 
events registered. For this purpose, we consider a total 
of 8 oscillation scenarios, which include the normal and 
inverted hierarchy. This demonstrates the range of possi- 
ble fluxes and energies that can be expected for different 
flavors. These estimates of the luminosity and energies 
of neutrinos from the disks of BH-NS mergers have a 
number of applications. 

Based in our results for the neutrino surface temper- 
atures, we estimate the number of events that would be 
registered in several detectors, both proposed and in op- 
eration if a BH-NS were to occur in the galaxy. While 
such events are expected to be rare in comparison with 
the standard core collapse supernova, it is still important 
to understand the potential signal. Our estimates can be 
used to rule out this object as the origin of a future neu- 
trino signal, as well as to guide investigations of neutrino 
detectors for the distant future. For these estimates we 
have included corrections due to the BH gravitational 
field such as energy redshift and bending of the neutrino 
trajectories. In general we find that there will be more 
counts coming from an AD-BH that is formed in a neu- 
tron star - black hole merger than counts registered from 
a PNS. For example in SK we will register 1000 more 
counts if a NS-BH merger happens than if a supernova 
event occurs. However, the neutrino signal from AD-BH 
would last only 0.15 s compared to the 10 s estimate for 
a PNS. We have then a signal that is 100 times more lu- 
minous and 100 times shorter than a PNS. We warn that 
our estimated counts are sensitive to the time the sig- 
nal lasts. And our estimate for the duration depends on 
the efficiency of converting gravitational energy into neu- 
trino energy. Therefore, the total counts could change as 
much as a factor of two. Nevertheless, our results for neu- 
trino energies, luminosity and fluxes are independent of 
this uncertainty. The amount of AD-BH neutrinos that 
could be detected in large scale facilities and the possible 
distinction of this signal from a PNS opens the question: 
Is there a way to distinguish the neutrinos from a BH-NS 
merger and a supernova? The timing of the signal and 
the energies of the neutrinos will be different, so this is 
an important clue. Therefore, it is not likely that we will 
mistake neutrinos from a merger event as those from a 
core collapse supernova. How could we guarantee a de- 
tection of BH-NS neutrinos? So far with a large scale 
facility like UNO we could detect neutrinos from AD-BH 
as far 4 Mpc, which covers galaxies in the local group such 
as Andromeda. To detect a BH-NS merger in neutrinos 
and compete with a gravitational wave detector such as 
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Advance LIGO would take a three orders of magnitude 
larger detector. 
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